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Abstract 



We study the scattering of a single parton state with a multi-parton state to derive the complete 
results of perturbative coefficient functions at leading order, which appear in the collinear factoriza- 
tion for Single transverse- Spin Asynimetry(SSA) in Drell-Yan processes with a transversely polarized 
hadron in the initial state. We find that the factorization formula of SSA contains hard-pole-, soft- 
quark-pole- and soft-gluon-pole contributions. It is interesting to note that the leading order pertur- 
bative coefficient functions of soft-quark-pole- and soft-gluon-polc contributions are extracted from 
parton scattering amplitudes at one-loop, while the functions of hard-pole contributions are extracted 
from the tree level amplitudes at tree-level. Our method to derive the factorization of SSA is different 
than the existing one in literature. A comparison of our results with those obtained by other method 
is made. 



1. Introduction 



In scattering processes with a transversely polarized hadron in the initial state. Single transverse- 
Spin Asymmetry(SSA) relative to the spin direction can be nonzero. SSA has been observed in various 
experiments, a review about the phenomenologies can be found in [1]. Theoretically, SSA can be predicted 
with the concept of QCD factorization, if large momentum transfers are large. In the factorization of 
SSA nonperturbative effects of the transversely polarized hadron are factorized into matrix elements of 
\ the hadron. Therefore, it will provide a new way to study the inner-structure of hadron by studying SSA. 
In this work we study the collinear factorization of SSA in Drell-Yan processes. 

From general principles SSA can be generated if the strong interaction changes the helicities of 
hadrons in a scattering and the scattering amplitude has an absorptive part. In the scattering involving 
a transversely polarized heavy quark, the helicity of which is not conserved in QCD because of the heavy 
mass. The related SSA can be calculated with perturbative theory of QCD, e.g., in [21 [3]. For light 
hadrons in high energy processes, one can neglect the mass of light quarks. The helicity of a massless 
quark is conserved in QCD. But this does not mean that the helicity of a light hadron is conserved in 
QCD, because a light hadron is a bound state of light quarks and gluons and the helicity of a light hadron 
is not only a sum of helicities of light quarks. 

The collinear factorization for describing SSA has been proposed in [31[5]. With the collinear factor- 
ization SSA in various processes has been studied in El [H O [101 ttH ttSJ [13] . In such a factorization, the 
nonperturbative effects of the transversely polarized hadron are factorized into twist-3 matrix elements, 
or called ETQS matrix elements. Taking Drell-Yan processes as an example, SSA is factorized as a 
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convolution of three parts: The first part is the standard parton distribution function of the unpolarized 
hadron defined with twist-2 operators. The second part consists of matrix elements of the polarized 
hadron defined with twist-3 operators. The third part consists of perturbative coefficient functions de- 
scribing the hard scattering of partons. If the factorization can be proven, the coefficient functions are 
free from any soft divergence like collinear- and I.R. divergence. In this approach the effects of helicity-fiip 
are parameterized with twist-3 matrix elements, while the absorptive part of the scattering amplitude is 
generated in the hard scattering of partons. 
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Figure 1: The cut diagram for the differential cross section of hadron scattering where the lower hadron 
is transversely polarized. The broken line is the cut. The black boxes represent parton density matrices 
of corresponding hadrons, the gray box is the forward scattering amplitude of partons. 

A widely used method to derive the factorization of of SSA at the leading order of ag is the diagram 
expansion at hadron level. All existing results are derived with the method except those in |14 ^ ll5 t[T6lll7j . 
This method has been also used for analyzing higher- twist effects, e.g., in |18tll9j. We take the Drell-Yan 
process hA + hs I'^i' -|- X as an example to illustrate the method. In the process Ha is transversely 
polarized with the spin vector s_i. The spin-dependent part of the differential cross section can be 
given by Fig.l. In Fig.l., the gray box represents Feynman diagrams for various contributions of forward 
parton-scattering with the cut. The lower black box represents the density matrix of the polarized hadron 
defined with quark- and gluon fields, and the upper black box represents the quark density matrix of 
the unpolarized hadron. The three parts are connected with parton lines. A collinear expansion of the 
momenta carried by the parton lines is performed to pick up the leading power contributions. For the 
partons from Ha the momentum is expanded around Pa, while the expansion for the partons from hs is 
around Pb- After the expansion, one obtains an approximated form of density matrices parameterized 
with various nonperturbative functions, i.e., parton distributions functions and twist-3 matrix elements, 
and the perturbative coefficient functions of the factorization. It is interesting to note that SSA in the 
factorization contains not only the so-called hard-pole contributions in which all three patrons from Ha 
carry nonzero momentum fraction, but also the so-called soft-pole contributions in which one of the three 
partons can have zero momentum fraction. 

It should be noted that QCD factorizations, if they are proven, are general properties of QCD. These 
factorizations hold not only at hadron level but also when one replaces the hadron states with parton 
states. The perturbatively calculable parts in factorizations, i.e., the perturbative coefficient functions, 
do not depend on hadrons and are completely determined by the scattering of partons. To derive the 
factorization of SSA, one can replace hadrons with parton states and calculate SSA perturbatively. The 
relevant twist-2 and twist-3 matrix elements can also be calculated with the parton states. In general the 
obtained results will contain soft divergences which usually appear beyond leading order. By writing SSA 
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as a convolution of these matrix elements and perturbative coefficient functions, one can determine the 
functions. In this work we will take this approach to derive at leading order all perturbative coefficient 
functions appearing in the collinear factorization of SSA in Drell-Yan processes. 

We notice that the approach taken here has been used to study factorizations only involving twist-2 
operators. Applying the approach for SSA, i.e., factorizations with twist-3 operators, will provide an 
alternative way to derive the factorization or to calculate perturbative coefficient functions. This will 
also give an independent verification of results derived with other approaches. It is not the intention 
here to prove the factorization or that the perturbative coefficient functions are free from soft divergences 
at any order. This is beyond the scope of the present work. However, as we will see, we already have 
at leading order some perturbative coefficient functions obtained by subtraction of collinear divergences 
with twist-3 matrix elements, in contrast to the case only with twist-2 operators. 

For the factorization only involve twist-2 operators, e.g., for the unpolarized differential cross-section 
one can simply replace each hadron state with a single parton state to derive the factorization. But 
for SSA, because of the helicity conservation of QCD it is not possible to obtain nonzero SSA and the 
relevant twist-3 matrix elements by replacing the transversely polarized hadron with a single quark state . 
But, one can construct multi-parton states to replace the polarized hadron. With the multi-parton states 
SSA and relevant twist-3 matrix elements are nonzero, because the helicity-ffip effects can be generated 
through correlations between these partons. 

In \14:\ [15] we have used multi-parton states to study the factorization. We have found that with 
tree-level results of SSA and twist-3 matrix elements there are only the hard-pole contributions. Later, 
in [15] it has been realized that there is a special class of one-loop contributions to SSA which can not be 
factorized as one-loop corrections to the hard-pole contributions at tree-level. These one-loop contribu- 
tions can only be factorized with some special twist-3 matrix elements at one-loop. These contributions 
are just the so-called soft-pole contributions. Their perturbative functions, although extracted from par- 
ton scattering amplitudes at one- loop, are at the same order as the hard-pole contributions derived from 
tree- level amplitudes. 

In this work we will use multi-parton states to derive all contributions in the factorization formula for 
SSA in Drell-Yan processes. They are Hard-Pole(HP) contributions, Soft-Quark-Pole(SQP) contributions 
and Soft-Gluon-Pole(SGP) contributions. There are two types of SGP contributions. One is from the case 
as given in Fig.l, where the gluon from the polarized hadron has zero momentum fraction. Another is 
from the case where three gluons are from the polarized hadron and one of them carries zero momentum. 
The twist-3 matrix elements for the three gluon case have been defined in {20) . 

It should be mentioned that besides the collinear factorization, there is another factorization for SSA 
in limited regions of kinematics. If the transverse momentum of the lepton pair is small, one can use 
the Transverse-Momentum-Dependent(TMD) factorization. The TMD factorization for unpolarized cases 
have been studied in [231 1211 [251 ESI [27]. For SSA nonperturbative effects of the polarized hadron are 
factorized into Sivers function[2l]. The properties of Sivers function and SSA with it have been studied 
extensively [28l|29l[30l[3ll[32l[33l[3ll[35l[36l[33. In[lll[l6l[I7]we have also examined the TMD 
factorization of SSA with parton states and found an agreement with existing results. 

Our work is organized as the following: In Sect. 2 we give our notations for Drell-Yan processes and 
the definitions of relevant twist-3 matrix elements. In Sect. 3 we introduce our multi-parton states. With 
these states one can define corresponding spin-density matrices in helicity space. The non-diagonal parts 
of the matrices are relevant for calculating SSA and twist-3 matrix elements. In Sect. 4 we study SSA in 
the scattering of multi-parton state at tree-level. With tree-level results we can derive HP contributions. 
In Sect. 5 and Sect. 6 we consider SSA at one-loop level. We find a special class of one-loop contributions 
which give the SQP- and SGP contributions. Sect. 7 is our summary. 
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2. SSA in Drell-Yan Processes and Definitions of Twist-3 Matrix Elements 

We will use the light-cone coordinate system, in which a vector a'^ is expressed as = {a'^ ,a~ ,d_i) = 
((a° + a^)/\/2, (a° - a^)/V2,a^,a'^) and = (a^)^ + (a^)^. Other notations are: 

with the light-cone vectors / and n defined as /'^ = (1,0,0,0) and = (0,1,0,0), respectively. We 
consider the Drell-Yan process: 

hA{PA, s) + hBiPs) ^-/*{q)+X +£+ + X, (2) 

where Ha is a spin- 1/2 hadron with the spin- vector s. We take a light-cone coordinate system in which 
the momenta and the spin are : 

(P+, 0,0,0), P^pa (0,P^,0,0), sf' = {0,0,s±). (3) 

The mass of hadrons are neglected. The spin of hs is averaged. The invariant mass of the observed 
lepton pair is = q^. We are interested in the spin-dependent part of the differential cross section, 
which can be written as: 

do' da / ^ N SyrOp q,« 

in which S = 2P^P^ . We parameterize the momentum of the lepton pair as: 

q^' = {xP^,yPB,q±). (5) 
The structure function Wr(x,y,gx) is related to the spin dependent part of the hadronic tensor 

^'^^ = E / -^e^'-^hA{PA,s±),hB{PB)\mrQm\X){X\qix)rq{x)\^^^ (6) 



X 

by: 



(-5/- + ^) = ef ^±a9±/3WT + • • • , (7) 

where • • • stand for spin- independent part. 

For large the structure function Wt can be factorized in the form of a convolution of perturbative 
functions with the standard parton distribution functions of and twost-3 matrix elements of Ha- The 
definitions of standard parton distribution functions with twist-2 operators can be found in literature. 
Here we discuss the definitions of twist-3 matrix elements of the transversely polarized hadron. The 
quark-gluon twist-3 matrix elements have been introduced in jU [5] firstly. We take a variant form and 
define them in the light-cone gauge n ■ G = : 

T,^{xi,X2) = .^,|^^e-^^(---)^"-^-^^"(P,.-jV;(yin) 
•7+ (G'+'^(2/2n) + i75G+^(y2n)) m\P,SL), 

r,_(xi,x2) = .x,/^e-^^^(---)^"-^^-'^^"(p,.-'Jv;(yin) 

•7+ [G+''iy2n) - i75G+^(y2n)) ^(0)|P, s^), (8) 
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with G+f" 



G^^. In other gauges gauge hnks along the direction n should be added to make the 



definitions gauge- invariant. One can also use the projection 7"*" or 7 75 to defined twist-3 matrix element 
Tqpixi, X2) and TgA.FC^^i, X2), respectively, as in [4J. The relations between these twist-3 matrix elements 
are: 

Tg+{xi,X2) = TqF{xi,X2) + TgA,F{xi, X2) , Tq^{xi,X2) = TqF{xi,X2) - TqA,F{xi, X2) ■ (9) 

One can show that the function TqF{xi,X2) is symmetric in xi and X2 and Tq/\^F{xi,X2) is anti-symmetric 
in xi and X2- 

The twist-3 matrix elements Tq±{xi,X2) with xi^2 > describe the correlation of those partons from 
Ha , which enter a hard scattering, e.g., the gray part of Fig.l. In the hard scattering, the initial quark 
carries the momentum faction X2, the gluon carries the momentum fraction xi — X2 and the final quark 
carries the momentum fraction xi. If the gluon momentum fraction is xi — X2 = 0, the corresponding 
hard scattering introduces the SGP contribution to SSA. If a quark carry zero momentum, i.e., xi = 
or X2 = 0, the corresponding hard scattering introduces the SQP contribution to SSA. It is clear that 
the SGP contributions are related to Tg+(x,x) with Tg+(x,x) = Tg_(x,x) = TqF{x,x), while the SQP 
contributions are related to rg-|-(0, x) or Tq±{x, 0). There are contributions with nonzero xi^2 and xi 7^ X2. 
These contributions are HP contributions. For the case xi < or X2 < the corresponding quark fields 
in the definition represent antiquarks. 

Instead of two quarks combined one gluon entering the hard scattering, there can be three gluons 
entering the hard scattering [20 j. The corresponding contributions can be factorized with matrix elements 
defined with twist-3 gluonic operators. In this case, as we will shown, there is a leading contribution of 
Os in the factorization of SSA. The contribution is a SGP contribution in which one of the three gluons 
carries zero momentum fraction. In general there are two types of twist-3 gluonic operators, distinguished 
by the color structure. One can define them in the gauge n • G = 0: 



0"'^^(xi,X2) 
7V°^^(X1,X2) 



9s ,bca f <iyi<iy2 iy-^xiP+ -iy2(x2-xi)P+ 
P+ J 47r 

{P, 5x|G^'+^(yin)G^'+^(y2n)G"'+"(0)|P, s^) 

■ 9s fbca f ^^;y}dy2_ iyj^xiP+ -iy2ix2-xi)P+ 

'P+^ J Att 

{P, sx|G^'+^(yin)G^'+^(y2n)G"'+°(0)|P, s^) 



(10) 



There are two scalar functions can be defined for each type of color structure in the case of xi = X2 = x 
for SGP contributions: 



0"^^(x,x) =gfPxGdi{x) + 

iV"/37(^^^) =gfpxGfi{x) + 



'9Ts^ + 9l's'^ 



xGd2{x), 
xG/2(x). 



In general the function Wt in the collinear factorization can be divided into four parts: 



Wt = Wt 




+ Wt 


+ Wt 






HP 


SQP 


SGPF 


SGPG 



(11) 



(12) 



For SGP contributions one can have two different contributions factorized with the quark- gluon- or purely 
gluonic twist-3 matrix elements, denoted by the subscriber SGPF and SGPG respectively. Each of the 
four parts can be expressed as convolutions of parton distributions, twist-3 matrix elements discussed in 
the above, and perturbative coefficient functions. Details of the convolutions will be given in the following 
sections. 
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The goal of our work is to derive all perturbative functions at leading order of a^. For HP contributions 
at leading order we only need to calculate with parton states parton scattering amplitudes and twist- 
3 matrix elements at tree-level. At one-loop level, there are in general collinear divergences in Wt- 
As observed in {15], at one-loop there is a class of contributions, whose collinear divergences can not 
be subtracted by using one-loop results of twist-3 matrix elements in the factorized HP contributions 
derived with tree-level results. These contributions hence can not be taken as one- loop corrections to 
the perturbative coefficient functions in HP contributions. In fact, the collinear divergences can be 
subtracted by the so-called soft-pole twist-3 matrix elements in which one parton carries zero momentum 
fraction. This is the origin of the soft-pole contributions. The soft-pole matrix elements are zero at tree- 
level, but nonzero at one-loop. This results in that the perturbative coefficient functions of the soft-pole 
contributions are at the same order of those of the HP contributions derived from tree-level. 



3. Spin-Density Matrices and Multi-parton States 



We consider a system |A^[A]) with total spin 1/2. The system moves in the z-direction with the helicity 
A = ± and can be a superposition of various multi-parton states. We consider a forward scattering of 
the system through some operator O which do not change helicity of quarks. The operator can be those 
used to define twist-3 matrix elements, or the hadronic tensor. In the later, the forward scattering is 
with some additional particles which are unpolarized. The transition amplitude is given as: 



Mx,x, = {N[X2]\0\N[Xi]). 



(13) 



We use Ai^2 = ± to denote the helicity of the initial- and final state, respectively. The transition amplitude 
in the helicity space is 2 x 2 matrix and can be expanded as: 



-MA2A1 = a + b-a 



A2A1 



or it can also be described with a spin vector = {s ,s): 

M{s) = a + b-s , 



-1. 



(14) 



(15) 



Prom the above the transverse-spin dependent part is determined by i.e., the non-diagonal part in 
the helicity space. SSA appears if the non-diagonal part of the hadronic tensor in Eq.dG]) in the helicity 
space is non zero. 

Because of helicity-conservation of QCD the non-diagonal part of M in Eq. (ll4p is zero, if |iV[A]) is a 
single quark state. Instead of a single quark one can consider the following multi-parton state: 



= \q[\]) + ci\qg[\]) + C2\qqq[\]) + Cz\qgg[\]) + 



(16) 



where all partons move in the 2;-direction, the sum of helicities of partons is 1/2 or —1/2. We will give 
later the details about the momenta and color structure of these partons. The • • • stand for possible states 
with more than 3 partons. We do not need to consider the states with more than 3 partons. The reason 
is the following: The leading power contributions to SSA come from those parton scattering processes, in 
which only three patrons from the polarized hadrons are involved. We call these involved partons active 
partons. Certainly, there can be more than three partons as active partons. The resulted contributions 
are power-suppressed and may be factorized with operators whose twist is larger than 3. In the leading 
power contributions, the three active partons can be the combinations of quarks and gluons. They are 
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995) Q<i9 and ggg. Charge-conjugated combinations should also be included. In the case described by 
Fig.l. the active partons are qqg. It is clearly that with first three states in Eq. (jl6p one can have all 
combinations by taking some partons as spectators. 

If we calculate the non-diagonal part of M in Eq. (I14p with the multi-parton state in Eq. ()16p , one will 
find with the helicity conservation of quarks nonzero contributions only from the interference between 
different states in the right hand side of Eq. (fT6|) . If we replace Ha with the above state and Hb with a 
single unpolarized parton in Eq.([6]), we will also get nonzero result for the spin-dependent part of W^'^ 
or for Wt. Similarly, the defined twist-3 matrix elements are also nonzero with the state |A^[A]). These 
nonzero results allow us to study the factorization of SSA. 

In |161 [T7] factorizations of SSA in Drell-Yan processes have been studied with the first two terms in 
Eq. (116p in the kinematical region q^/Q"^ <C 1. In this case, all partons are active. Non of them can be a 
spectator parton. But for interferences between other states, some partons can be spectators, because we 
only need to consider those interferences with three active partons. The existence of possible spectators 
only affects overall factors of interested quantities like Wt and twist-3 matrix elements, it has no effect 
on the derivation of perturbative coefficient functions. Hence, for our purpose we only need to consider 
those matrix elements {ab\0\c) or {c\0\ab), where a, h and c are quarks or gluons. These matrix elements 
can be obtained from the interferences between different states in Eq. (jl6p by taking out some partons as 
spectators. We will illustrate this in the following. 

For the interference between the q- and the g(7-component, non of partons can be a spectator. We 
define the state \q[X\) and the state |<?5[A]) as: 

\q{p, \)) = blip, Ag)|0), \q{pi,XMP2, Ap)) = Tl,/.^{pi, X,)ai(P2, A<;)|0), 

p^ = (p+, 0,0,0), pt = XQp^, p2 = (1 - xoK = xop'^, (17) 

where h\ is the quark creation operator with i as the color index, a\ is the gluon creation operator with 
a as the color index. \q{\g) is the helicity of the quark(gluon). To simplify the notations we will write 
p~^ = Pj" and p~ = Pj^. It is straightforward to obtain the non-diagonal part as 



{q{p, +)\0\q{pi, +)g{k, -)) + {q{p^, -)g{k, +)\0\q{p, -)) 
-)\0\q{pu -)9{k, +)) + {q{pu +)g{k, -)\0\q{p, +)) 



(18) 



In the above we use the index qg to denote this type of interference contribution. We also introduce a 
coefficient C^^ in the non-diagonal part of the spin-density matrix. Contributions of this type to twist-3 
matrix elements and Wt will be proportional to C^, and will be called g(7-contributions. The derived 
perturbative function will not depend on C^^ . 

For the contribution from the ggg-state we note that the interference with the single quark state is 
zero if the total helicity is changed. The interference with the qgg-state does not need to be considered 
because at least four partons must be active and one quark is a spectator. Hence, we only need to consider 
the interference with the g^-state, in which one quark is a spectator and three partons are active. In 
this case, the forward scattering is participated by a gluon and a gg-pair. In order to have AA = ±1, 
the total helicity A of the gg-pair must be zero. There can be two states with A = for the gg-pair. We 
denote the two states as : 

I iQQ)±) = Tl,^ +)<(P2, -) ± bl[pu -)dl{P2, +)] |0), (19) 
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and the single gluon state as \g{Xg)). The gluon carries the same color index a and the momentum p. 
Xg is the helicity. is the create operator for the antiquark. With these states one can construct the 
non-diagonal part of the spin-density matrix j\4^'^^^ as: 



M 



M 



m 

-+ 



cf 
cf 



((gg-)+|0|5(-)> + (9(+)|C?|(QQ-)+) 
((gg)+|0|5(+)> + (5(-)|0|('Z^-)+) 



m)-\0\9{^)) - l,g{-)\0\{qq)_ 



(20) 



In the above we have introduced two coefficients to distinguish the contributions from the two qq 
states. We note that there is a sign difference for the terms with Ci'' between the first and the second 
equation. This difference can be easily found by requiring that the state | {qq)_) becomes a spin-1/2 
system by adding a quark. Again, we will call all contributions from this matrix as gg-contributions. 
They are linear in the two coefficients Cj^ . The derived perturbative functions will not depend on C'^ . 



For the contribution from the g^^-state in Eq. (jl6p . only the interference with the g^r-state and with 
the qgg-state need to be considered here, where one quark can be taken as a spectator. In this case, one 
has the forward scattering as gg ^ g or g ^ gg. The color of the two gluon state must be the same 
as the color of the one gluon state. The total helicity A of the two gluons must be zero. There are two 
states with A = for a given color structure. We denote 

l(55)±) =ir''=[4(pi,+)4(P2,-)±4(pi,-)at(p2,+)] |0). (21) 

With these states one can construct the non-diagonal element of the spin-density matrix Jli^ds^^ as: 



M 



M 



(ggF) 



{ggF) 



■99 



■99 



((55)+|0|5(-)) + (9(+)|0|(59)- 
((55)+|0|ff(+)) + (5(-)|0|(99)- 



■gg 



■gg 



{{gg)_\0\g{-))-{g{+)\0\{gg)^ 
{{gg)_\0\g{+))-{g{-)\0\{gg)^ 



.(22) 



In the above we introduce two coefficients J-j^ to distinguish the contributions from the two states in 
Eq. (l2ip . Another spin-density matrix Jli^99^^ can be constructed in this case by replacing i/"'"^ in Eq. ()2ip 



with d"'^'^, and with V'^ in Eq.(|22p. We will call all contributions from these two spin-density matrices 



as ^'^'-contributions. They are linear in the four coefficients J-^ 



and . 



With the constructed spin-density matrices in the above one can calculate twist-3 matrix elements 
and the structure functions Wt by taking correspond operator O. It is straightforward to obtain the 
twist-3 matrix elements Tgj- at tree-level. The results from the qfj-contributions are 

rg+{xi,X2) = Ci^ngsV^^iN^ -l){x2-xi)5{l-xi)6ix2-xo) + Oigl), 



Tg+{Xl,X2) 
Tg^{Xl,X2) 



C^Hgs 

-C'?%5.\/2^(iVc - l)(a;2 - xi)5{l - X2)5{xi - xq) + 0{g\ 



(23) 



The results from the gg-contributions are 



Tq+{XI,X2) 



TTgs{N^ - l)\/2xoxo 



{cf 



C_) 6{xi + xo)S{x2 - xo) 



Cf + C'L')5{x2 + xo)6{xi-xo] 



0{g% 



Tq_{xi,X2) = irgsiN^ - l)\/2xoXo 



Cl' + CI') 6{xi + xo)6{x2 - xo) 



+ 



(Cf-Cl^) 5ix2 + xo)6ixi-xo] 



+ 0{gl). 



(24) 
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With Wt calculated in the next section at leading order, one can find the factorized form of Wt in terms 
of Tq± with the above results. This is for HP contributions. In Sect. 5 and 6 we will also give the results 
of Tq±{x,x), Tq±{Q,x) and these gluonic twist-3 matrix elements. These results are at order of and 
will be used to factorize the soft-pole contributions. 



4. Hard-Pole Contributions 



As discussed in the last section, we replace the polarized hadron with the multi-parton state in 
Eq.([T6]) to calculate the non-diagonal part of the constructed spin-density matrices for Wt- We replace 
the unpolarized hadron hs with single-parton states. In this section we will work at tree-level. 





R 
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Figure 2: The diagrams for the amplitude q + {q + G) — t- ^* + X — t- g + g at tree-level. The black 
dots denote the insertion of electromagnetic current operator. Broken lines represent the cut. For the 
propagators with a short bar only the absorptive part of the propagator is taken into account. 

We first consider the gf^-contributions. If we replace hs with an antiquark q with the momentum 
= (0,p~,0,0), the leading order contribution to Wt comes from diagrams in Fig. 2. The complex 
conjugated diagrams should be included in order to obtain the non-diagonal part of the spin-density 
matrix J\A^'^^^ given in Eq. ()18p . In the diagrams of Fig. 2. the broken line divides each diagrams into 
a left- and right part. Each part represents a scattering amplitude. The short bar cutting a quark 
propagator is in fact a physical cut of the amplitude represented by the left part. It means that only the 
absorptive part of the cutting propagator is taken into account: 



Abs 



il ■ kg 
/Cg + ie 



(25) 



It is straightforward to calculate these diagrams and we obtain: 
Wt 



Figm 



^,g^2 9sa^V^Nlj.^^_ _ y-^^^^^^^ -x){l-y)- ql){Nl + y - I) 



1 



1-y 



4vr ql iV2 
y"^ + XQ- |Ag|(y2 - xo) 



(26) 



with s = 2p^p~ . Cq is the electric charge of the quark q in unit e. The 5- function of 1 — x — yxQ is from 
the cutting quark propagator. The terms with \Xq\ = 1 are quark-spin dependent, because the external 
quark lines are extracted with Xq^^'y ■ p. 
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With the gg-contributions of Tqj- one can write the above Wt into a factorized form. The terms with 
|Ag| = 1 should be factorized with Tq+ — Tg^ or Tq/\^p, because 7"'' 75 is used to define them. The other 
terms should be factorized with (Tq+ + Tg_) or Tqp. With Tq± in Eq. (j23p we have the factorized form: 



Figm 



Tr^Ncqj_ Jx yi Jy y2 

+'Hq-{^l,^2)Tq-{yi,XB 



'Hq+iCl,C2)Tq+{yi,XB) 



(27) 



with 



6 



2iVc6(i-ei)(i-6) 

X y q- 

— , 6 = — , XB = 

yi y2 ^q-p 



(6 + 6 

,2 



1), 



N^ + ^2-i 



2Afe(l-6)(l-e2 



^2 



yiy2s. 



(28) 



The function fq{y2) is the antiquark distribution function of hs- For Hb = q, we have fq{y2) = <5(1 — 
2/2) + C'(as)- It is noted that the derived perturbative coefficient functions do not depend on C^^ . One 
can also replace Hb with a quark. In this case the results can be obtained by reversing the directions of 
quark lines in Fig. 2. They can be obtained from the above results through charge-conjugation. We will 
give them at the end of this section by combining all parton flavors. 
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'(xaaaoaoa 
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■000000000 
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(d) 




Figure 3: The diagrams for the amplitude G + [g + G] — )• 7* + X — )• G + 5 at tree- level 

If we replace the unpolarized hadron Hb with a gluon carrying the momentum p, the leading contri- 
butions to Wt comes from Fig. 3. The calculation of these diagrams is similar to the calculation of Fig. 2. 
We have the sum of Fig. 3: 



Wt 



Figm 



2 aso-s 



^2x^6{s{l - x){l -y) - q\)5{x - yxo) 



l + jy- 1)N^, 



xo{l-yf + y^ + \Xq\{xo{l-yf-y^ 



(29) 



This result can be factorized in the following form: 



Wt 



/ / fgiy2)6{s{l - - C2) - q±) 

Jx yi Jy y2 



Figm 



yi Jy y2 
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1 + (6 - i)iv<' 



2(iV2_i)(i_^^)^2 
l + fe-l)A^c' ,2 
2(iV2-l)(l-6) 



(i-6)'(i-6-e2), 

S2- 



where fg{y2) is the gluon distribution function. For hs = g we have /g(y) = S{1 — y) + 0[as 






' 1 


• — > — 






mm 


• 1 ^ 


[c) 



(d) 



(30) 



Figure 4: The diagrams for the amphtude q+[q + q\^^* + X^q + GaX tree-level. 

By replacing hs with an antiquark q{p), the gg-contributions for Wt are also at leading order. They 
are given by the diagrams in Fig. 4. In the first four diagrams the anti-quark q{p) in the initial single 
parton state must have the same flavor as the quark in the multi-parton state, while in the last four 
diagrams q{p) can have different flavor. The results are: 



2 asOis 



^^^^ 



(C- + - 2^o)) <^(s(l -x){l-y)- ql) 



'5{x - yxo) 



■y 



|l4l4]f4|4K4|l 



1-y 

2 9sOLs 



-iV2 - 1 



{'^-yf^x-yx^) , 

(Ci'? + (1 - 2xo)) <5(s(l - x)(l - y) - ql) 



b{x- yxo) + 5{x- yxo) 



(y2 - 2y + 2 



(31) 



With the tree-level results of the gg-contributions for the twist-3 matrix elements we can derive the 
following f actor ized form: 



+ 



eqOs ^ dyi ^ dy2 2n 
T^Ncq\_ Jx yi Jy y2 

C2)Tg+{-yi, Xb) + 'Hgq-{£.l,(.2)Tg-{-yi, Xb) 
'^qq+{iui2)Tq+{-XB,yi) + Tigq-i^l, ^2)Tq-{-XB , IJl) 
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dyi dy2 



Tr'^Ncql Jx yi Jy y2 

+-^5^0(6,^2) 



fg,{y^)6isil-^l)il-C2)-ql) 
Tq+{-yi,XB) - Tq^{-XB,yi) 

Tg+{-XB,yi) - Tq^{-yi,XB) 



(32) 



with yi = yi — xb and the perturbative functions: 

1-6-6 



^99+ (6 16) 

^gg+(6)6) 

"^990(616) 



2iVc6(i-6)' 
(1-6)^(1-6 



) 'Wqg-(6>6 



1-6 



2iVc6 
-26 + 2 

2ei 



2A^c6(l-6)' 

^g-g-(6,6) = (i-6-6) 



(1-6-6), ^5,0(6,6) 



(1 - 6)^ 

2iVc6 ' 
il -26 + 2 

2ei 



(1-6 



(33) 



fq'{y2) is the antiquark distribution function for the flavor which does not need to be the same as the 
flavor of quarks used to calculate the twist-3 matrix element Tqj-. 

The studied contributions plus charge-conjugated processes give the all leading HP contributions for 
SSA. All perturbative coefficient functions are at order Og. Combining all possible apron flavors we obtain 
the factorized HP contributions as: 



Wt 



a. 



HP 



1 dyi dy 



iT^Ncq\ Jx yi Jy y2 



■5(^(l-6)(l-6)-'zi) 



'Wg+(6, 6) (^lfq{Vi)Tq+{yi,XB) + ^g-(6, 6) Y ^^g fqiy2)Tq-{yi, Xb) 

[q] [q] 

+'Wg5+(6,6) XI ^g-^9(2/2)Tq+(-yi, xb) + -H<;5_(6, 6) XI ^-b) 
[q] [g] 

+^gg+(6, 6) X ^9/9(y2)^g+(-a;B, yi) + ^gg-(6, 6) X ^?./'5(j'2)T5_(-a;B, yi) 



+^gg0(6,6) X ^Ifg'iy^) 
[<?,<?'] 

+^gg0(6,6) X ^Ifg'iy^) 



Tq+{-yi,XB) - Tq^i-XB,yi) 
Tq^{-XB,yi) - Tq^{-yi + xb,xb) 



+^g+(6, 6) X (^lfa(y2)Tq+{yi,XB) + ^g-(6, 6) X eg/(;(y2)rg_(yi, xb) 

[g] Ig] 
where the notation for summing over flavors is defined as: 



(34) 



q=u,d,s,- 



X eg/g'(y2)rg±(2;i,Z2) 



E 



fqiy2)Tq±{zi,Z2) - fq{y2)Tq^{- Z2, -Zl) 

fq'{y2)Tq±{zi,Z2) - f q' {y2)Tq^{- Z2 , - Zl) 



[9,-?'] 



q=u,d,s,---,q'=u,d,s, 



J2'4f9iy'2)Tg±i^i^^2) = E 6^/9(^2) 

[9] q=u,d,s,--- 



Tq±{zi, Z2) - Tq^{-Z2, -Zl) 



(35) 
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It is interesting to study the limit q^/Q"^ <^\hy using 



5(1-6) , 5(1-6) 



+ 



(1-6)+ (1-6)- 

In this hmit, the above contribution to Wt becomes: 



5(l-6)<5(l-6)ln^. 



(36) 



HP 





f6{l 




Jy y2 


1(1- 


ei)+[ 



6 Yl eg/g(y2)Tg+(yi, x) + elfqiy2)Tq-{yi,x) 

[q] [</] 



+5(1-6 
, 5(1-6) 



1 + 6^ 
(1-6)- 



1 + 



6-1 



25(1 -6) In 



H ^9/5(^2)^9+ (yi>yi) 
[9] 



+ 



iV,(iV2 _ 

■^(1-6) 
m 



'-{N^{1 - 6) - i)((i - 6)' + 6') E elfg{y2)Tq+{yi,yi] 

' \1\ 



(1 - i\)Y^\hiV'iy^q-^-V^ + X,X) - i\Y^qh'KV2)Tq+{.-y\ + X, x) 

[9] [9] 



i + o((?i/g2 



(37) 



It is noted that in the hmit SGP contributions appear. If we take the limit q^/Q"^ ^ 1 in the tree-level 
results for the parotnic Wt's in Eq. ()26p and Eq. ()29p instead of the factorized results in Eg. (1271) and 
Eq. (j30p . we will not obtain the SGP contributions. However, the SGP contributions can be derived by 
using parotnic Wts in the limit beyond tree-level |17j. 

The factorized results of Fig. 2 and Fig. 3 have been derived in [11] with the method of the diagram 
expansion mentioned in the Introduction. By rewriting the above results with partonic variables which 
are defined as: 



t 
u 



{viPa - qf = -s6(i - 6) - s = yms, 



{y2PB 



-s6(i-6)-gi, Q' = g' = 66S-(?i, 



(38) 



we find that our results agree with them in [llj . Recently, the results corresponding to the contributions 
from Fig. 4 have been derived with the method of the diagram expansion in [8j . Again our results in Eq. 
([32l) agree with those in [8]. 



5. Soft-Gluon-Pole Contributions 

The SGP contributions comes from the case when one gluon with zero momentum enters hard scat- 
tering. They may come from the qg-, qq- and the ^gi-contributions. The qg- and qq contributions 
are factorized with the quark-gluon correlator Tq^{x,x) = Tq-{x,x). Later we will show that the qq- 
contributions need not to be studied, because it is automatically included in the factorized form obtained 
from the g^-contributions. The ^^-contributions are factorized with the purely gluonic correlator defined 
in Eq. (112p . We will discuss these two types of contributions in this section separately. 

5.1. The g(jr-Contributions 

We have given the results of the g(7-contributions for Tq±{xi,X2) at tree-level in Eq. (j23p . At this 
order one simply has Tq±{x,x) = 0. However, beyond the tree-level, Tq±{x,x) can be nonzero. As found 
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Figure 5: The diagram for Tg±{x,x) at one-loop. The black dots denote the insertion of operators used 
to define Tq±{xi, X2) in Eq.(l8|). 



in I17j . at one- loop level there is only one diagram giving nonzero contribution to Tq±{x,x) in the 
light-cone- or Feynman gauge. The calculation of the diagram is straightforward. The contribution has 
an U.V.- and a collinear divergence. Both are regularized with the dimensional regularization as poles of 
e = 4 — d. After extracting the U.V. pole we have|161 117]: 



r,±(x,x,M) = -C^^N.iN^ - l)xo 



x) 



ecJ 47r/x^ 



(39) 



where the pole is the collinear divergence with the index c. /i is the renormaliation scale related to the 
U.V. pole, and fic is that related to the collinear pole. 

To find out the SGP contributions it is convenient to work with the light-cone gauge n ■ G = 0. We 
consider a special class of diagrams which represent a part of one- loop corrections to those given in Fig. 2. 
These diagrams are obtained from Fig. 2. by adding a gluon. They are given in Figl6j In the first four 
diagrams the gluon is emitted by the initial gluon and is absorbed by the final quark. In the last four 
diagrams the initial gluon goes across the cut represented by the broken line and emits a virtual gluon 
absorbed by the outgoing quark. 

The contributions from Figl6] contain a collinear divergence. In the first four diagrams, the divergence 
appears when the lowest gluon crossing the cut is collinear to the +-direction. In the last four diagrams it 
appears when the gluon emitted by the outgoing quark is collinear. Because the contributions from Figl6] 
are one-loop corrections to Fig. 2, one may expect that the collinearly divergent parts of the contributions 
can be re-produced in the factorized form of the contributions from Fig. 2. in Eq. (l27p . where one replaces 
Tq±{xi,X2) with the corresponding one-loop Tg-|-(xi, X2). As discussed in detail in [T7], this is not the 
case, because the color factor here does not match. Even if one neglects the color factor, the divergences 
still can not be re-produced. 

Analyzing the collinear divergences in the contributions of Figl6l one finds that the collinear diver- 
gences should be factorized with Tq+(x,x) = Tq_(x,x). Taking Figl6^ as an example, the added gluon is 
with momentum ki. If ki is collinear to the +-direction, i.e., k^ ~ (1,A^,A,A) with A <C 1, one can find 
that the gluon exchanged between the initial gluon and the initial antiquark is soft with the on-shell con- 
dition of the cut propagator. In fact, this gluon is a Glauber gluon with the momentum ~ (A^, A^, A, A). 
Comparing Figj6^ with Fig. 5, one can identify that the gluon crossing the cut in Fig. 5 corresponds to the 
collinear gluon with ki in Figj6^. If the collinear gluon is contained in Tq±, the Glauber gluon should be 
taken as the gluon entering hard scattering. Since it is a Glauber gluon with vanishing momentum, the 
divergent parts of Figj6]should be factorized with Tq^{x, x). This is the reason why the SGP contributions 
appear. 
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Figure 6: The diagrams for the amphtude q + [q + G] -^^*+X^q + q for SGP contributions. The 
black dots represent the insertion of electromagnetic current operator. 



Performing the same analysis for FiglBb, FiglGt and Figl6]i in the case that the gluon crossing the cut 
is collinear, one will find that the gluon exchanged between the initial antiquark and the initial gluon is 
a Glauber gluon. For the last four diagrams the gluon emitted by the outgoing antiquark in the right 
part is a Glauber gluon, if the gluon emitted by the outgoing quark is collinear. Therefore, the collinear 
divergences in these diagrams are related to the Glauber gluon. It should be noted that only the diagrams 
in FiglU contain such a collinear divergence related to a Glauber gluon. 

Before giving the results, the following facts should be pointed out. In Feynman gauge, one has to 
consider more diagrams which contain the collinear divergence, e.g., instead of that the collinear gluon 
is attached to the initial gluon in the left part of FiglU the gluon can also be attached to the initial 
antiquark. Such diagrams are finite in the light-cone gauge, at least for most cases studied here with 
an exception which will be discussed in Sect. 6. In the following we will work in the light-cone gauge 
n-G = 0. 

The contributions of Figl6]contain an integration of a loop-momentum. It is easy to find the collinearly 
divergent part of the contributions by expanding the integrand in A, where the collinear gluon has the 
momentum ~ (1, A^, A, A). We find the collinearly divergent part of the contributions from FiglUas: 



27r2 2Nr 



Xo 



— Ixxq — x\y 
(xo - - y)s 



x^ + Xgy^ 
(xo - x)(l - y) 



6{u) = 6{s{xo - x){l -y)- q\). 



(40) 



In the above the pole in Cc = 4 — d represents the collinear divergence. The (5-function from the on-shell 
condition of the intermediate gluon exchanged between quarks also depends on the loop momentum and 
needs to be expanded in A. This results in the terms with the derivative of the (5-function. The last four 
diagrams do not contain terms with the derivative of the 5-function. With the result of Tq^{x,x) from 
the (^^-contribution in Eq. (|39p we can derive the factorized form: 



Wt 



e^as 



FigM 



IT 



2 AT, 



cqi 



Jx yi Jv 1/2 
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5Gg(6,C2) yi 



dyi 



Nc{l - 6) ' 



2gi(i-6)'-gi-eie2(2-6: 

iVe(l-ei)(l-6) 



(41) 



We note that the perturbative coefficeint function here is at the same order of as as those of HP contri- 
butions because Tg+{yi,yi) is at the order of gsCtg. 
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Figure 7: The diagrams for the amphtude g + [q + G] -^'y* + X^g + q for SGP contributions. The 
black dots represent the insertion of electromagnetic current operator. 

If we replace Hb with a gluon, one obtains similar diagrams from the SGP contributions from the 
ggf-contributions. These diagrams are given in FigiTl The collinearly divergent part of the contributions 
belong to the SGP contributions. We have calculated the collinear divergences in these diagrams in 
the light-cone gauge and in Feynman gauge. The same results are obtained. This corresponds to the 
situation with Tq±{x,x) with Figj5l only the same one diagrams in the two gauges gives the result in 
Eq.dM])- From Figdwe have: 



Figm 



47r2 ^ 
, s{l - y)S{u) 



^^0 



s(l - + 2xxoiy " 2) + 2x^(^2 _ 2y + 2)) 



6'{u) 



2 ^2 



{xo - x){-xy + 2x + 3xoy - 4xo) - xly{y^ + (1 - yf)^ ' 



(42) 



With the result of Tq^{x,x) from the g(7-contribution in Eq. (|39p we can derive the factorized form: 

'■1 dyi /"i dy2 



ir^Ncqi Jx yi Jy y2 



Ql) 



c tc c \ l ^Tq+{yi,yi) 

5g9(6,6) yi 



dyi 



+ SGg{(l,(2)Tq+{yi,yi) 



5G3(ei,6) 



iV2 



J (e? + 26(6 - 2) + 2(6' - 26 + 2)) , 
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5G<,(ei,6) 



rice _LfiC 0C2 , Of A 6(^2 + (1 - 6)^) 

34i?2 + 6^1 - + 3^2 - 4 



The factorized results have also been also derived with the method of diagram expansion in 



(43) 



p 

f 0(55530' 



Figure 8: (a). The possible SGP-contributions from the gg-contributions. (b). The diagram for Tq-|-(x, x) 
in the gauge n ■ G = from the (/^-contributions. See the discussion in text. 

In the case when hs is replaced by a gluon, one can have the SGP contribution from the qq- 
contributions. An typical diagram is given in FiglHl One can also obtain Tq±{x,x) from the qq- 
contributions at this order. The diagram for it is given by Figilb. It is easy to find that the SGP 
contribution is included in the factorized form in Eq. (j43p . 

Combining contributions of all flavors the SGP contributions can be factorized with the quark-gluon 
twist- 3 matrix element as: 



a. 



SGPF TT^Ncqi Jx yi Jy y2 



[1] 



[9l 



.(44) 



The above results agree with those in |11|,[8] derived with other method. Again, in the case of q^/Q"^ <^ 1 
this contribution takes a simplified form: 



a.. 



SGPF 7r2Ar2(g2 )2 J ,^ Jy 



dyi dy2 \, .2^xn c ^ 2^/ dTq+{yi,yi)\ 
Jx yi Jy y2 V ^ ^ V oyi J 



+ 



i^^(2e? - 3e? - 1) - #M^6(i + il) + m - 6)5(i - 6) in ^ 



(1-6)- 



(1-6)+ 



■Y,elfgiy2)Tq+{yuyi) + ^^2^(6' + (1 - 6)')<5(1 - 6) E e'^(y2)r,+ (yi, yi) 

M [q] 



+ ■■■ 

where • • • stand for contributions suppressed by q'^/Q'^- 
5.2. The (/(/-Contributions 



(45) 
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At the order we consider, there is no HP contribution from the (jr(7-contributions. But, it is possible 
that there are leading SGP contributions from Wt at one-loop level, similar to cases considered in the 
above. We consider first the gluonic twist-3 matrix elements in Eq. (|lip . These functions are zero at 
tree-level. 




(6) 



Figure 9: The diagrams for G(iiji{x) and Gd2j2ix) in the light-cone gauge. 



At one-loop level, the functions become nonzero. They receive nonzero contributions from the di- 
agrams given in Figl9] in the light-cone gauge. In Feynman gauge there are more diagrams. In this 
subsection we will work with the light-cone gauge. For the factorization studied below we only need to 
calculate Figj9^ and the corresponding diagrams for Wy. The contributions from Figl9]3 and the corre- 
sponding contributions to Wt can be obtained from the permutation of the two initial gluons. We will 
only give results from FigJO^ and the corresponding results of Wt- We obtain: 



Gdiix) 
Gd2{x) 
Gfi{x) 
Gf2{x) 



gsOis 



2y/2 



{N^, - A){Nl - l)5{x - xo) 



+ In 



eV 



4^/2 



{Nl-A){Nl-l)5{x-xo) 



— I + ln 



2^/2 

gsds 

4^/2 



A^2(iv2 _ i)5{x-xq) 
N^iNl - l)5{x 



+ In 



xo) 



2 



+ In 



/i, 

/2, 



(46) 



with the parameters di^2 and /1.2 related to J-j^ in Eq.(|22p and T)^^ as: 

di = (1 - xo) (xoPf + V3J>) , d2 = + 



(47) 



The corresponding contributions to Wt are given by diagrams in FigllOl The results for the color 
antisymmetric gluon state are: 



Wt 





5'{u) 







-/2 (4x2 + 4xxo(y - 2) + 4{y^ - 6y + 6) 



2/1 (2x2 ^ 2xxo(y - 2) + xliy - 2) 
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Figure 10: The diagrams for the amplitude q + {G + G) — t- 7* + X — t- (7 + G at one-loop for possible SGP 
contributions. 



6{u) 



4/i(2x + xo(y-2)) + 2/i 



x{5y - 4) + xo{3y^ - 7y + A) 
1-2/ 



s(xo - x)(l -y) - q\. 



(48) 



Replacing the color factor N^{N^ — 1) with —{N^ — 4:){N^ — 1) and /i^2 with di^2i respectively, one obtains 
Wt from FigllOl with the color structure of d"'^'^. With the results of the gluonic twist-3 matrix elements 
in Eq. (j46p we can derive the factorized form from the SGP contribution from FigllOl bv combining all 
flavors as: 



Jx yi Jv y2 



j=l,2 



ff \i dGfi{yi,yi) dGdi{yi,yi) 
fq{y2) { yi — —_ h yi 



+/g(y2) yi 



dGfi{yi,yi] 



dyi 



yi- 



dyi 

dG di{yi,yi 
dyi 



dyi 



+ E '5G^(el,6)Ee 

i=l,2 q 

+fq{y2) {Gfi{yi,yi) - Gdi{yi,yi)) 
with the pertubative functions: 

'Sgi{S,i,S,2) 
<Sgi{S,i,S,2) 

<SG2iCl,C2) 



fq{y2) {Gfi{yi,yi) + Gdi{yi,yi)) 



^^(2^2 + 26(6 -2) + (6 -2) 

1 - ?2 ^ 



1-^2 



(4^? + 46 (6 -2)+ el -66 + 6), 



(ee? + 4^1 (26 - 3) + 36' - 106 + 8) , 



(l26' + 66(36 - 4) + 76' - 206 + 14 
i - 42 ^ 



(49) 



(50) 
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Prom the above results we can derive the result in the limit q± ^ as: 



SGPG 



with: 



TT^Nciqiy Jx yi Jy y2 
^ i=l,2 q 

,fr \( dGfi{yi,yi) dGdi{yi,yi] 

+fq{y2) yi — — yv 



. dGfi{yi,yi) dGdi{yi,yi] 

fq{y2) ( yi — — + yi 



dyi 



dyi 



dyi 



dyi 



i=l,2 



fq{y2) {Gfi{yi,yi) + Gdi{yi,yi)) 



+fq{y2) {Gfi{yi,yi) - Gdi{yi,yi)) 



S^i = 2e?- 2^1 + 1, 5^2 = 4e?- 4^1 + 1, 

s±i = - {efi - 4fi + 1) , = - (i2e? - ee? + 1 

The above the SGP contributions are leading contributions in the limit. 



(51) 



(52) 



6. SQP-Contributions 




Figure 11: The diagrams for the twist-3 matrix elements with X2 = in the gauge n-G = 0. The first two 
diagrams are for xi > in g(7-contributions. The later two diagrams are for xi < in gg-contributions. 



Similarly to the twist-3 matrix elements for SGP contributions, the twist-3 matrix elements for SQP 
contributions are zero at tree-level, because one can not define a quark state with zero momentum. 
Beyond tree-level, they can be nonzero. In the light-cone gauge n ■ G = 0, one can find two possible 
diagrams at one-loop for the ^^-contributions and the gg-contributions. They are given in Figllll It is 
easy to find that Figlllb and Figllld will give zero contribution. We have for the ^^-contributions from 
FigJTTb as: 
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Figure 12: The diagrams in the n • G = gauge which give the soft fermion pole contributions to SSA. 



We have for the (/(/-contributions from Figlllb as: 



T,+ (a;,0) 



{cf-C''_^)g,as5[x + xo 



N,{N^ - 1) 



V / 4 Xo 



1^ I i^^v 



(54) 



It is noted that in the above x is negative. It implies that an antiquark with the momentum fraction —x 
enters a hard scattering. 

The SQP contributions from the gg-contributions to Wt are given by diagrams in Fig|12lin the gauge 
n-G = 0. Following the analysis similar to that of Figl6l one can see that the vertical quark line in the left 
part of diagrams carries the momentum kg at the order of ~ (A^, A^, A, A), if the gluon at the bottom 
crossing the cut is collinear, i.e., its momentum scales like (1,A^,A,A). Factorizing the collinear gluon 
into the corresponding twist-3 matrix elements, one can realize that in the left part of diagrams, there 
is a gluon combined with a soft quark entering the hard scattering. Therefore, the collinearly divergent 
contributions are SQP-contributions. 

It is straightforward to find the divergent contributions from Figil21 



^qq "^Idsal ^/2^5{s{xQ - x){l - y) - q\) 



levr^A^s 



91(1 - xq) 



(x(2x-3(xo-l)(2/-2)) 



1 - Xo 



+(1 - xo)(y' - 5y + 5) - \\\{x{y - 2) - (xq - l){y^ -?>y + 3)) 



(55) 



Again the quark-spin independent part should be factorized with the combination Tj^q{x^Q) +T_„y(x,0), 
and the contribution with |Aq| should be factorized with T+q(x, 0) — T^q{x, 0). With the results in Eq. ([53l) 
we have: 



Jx yi Jv y2 



Fig[T2] 
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^^2^((ei + 6)' + 4(i-ei-e: 
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We turn to the (/(/-contributions. The contributions are given by diagrams in FigJTSl in the light-cone 
gauge. We need to find the collinear divergences related to the collinear gluon crossing the cut in these 
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Figure 13: The diagrams in the gauge n ■ G 
one-loop for possible SFP contributions. 



for the amplitude g + (g + g) — )■ 7* + X — (7 + G at 



diagrams. But, a direct calculation of the collinear divergences in these diagrams will give wrong results. 
This is the exception mentioned in Sect. 5.1 before Eq. (j40|) . We will explain this with Fig|13b as an 
example. In this diagram, the collinear divergence appears when the gluon attached to the initial quark 
is collinear to the +- direction. Instead of attaching the collinear gluon to the initial quark, it can also 
attached to other places. There are two examples given by the diagrams in FigJHl 
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Figure 14: The diagrams obtained from FigllSb by changing the attachment of the collinear gluon. 

As discussed in Sect. 5.1., one may expect that these two diagrams in Fig|14ldo not have the discussed 
collinear divergence in the gauge n-G = 0. Because of the structure of the color factor, Fig|14b^ is always 
zero. But, through an explicit calculation one finds that Figll4b also contains the collinear divergence. 
Similarly to Fig|13b. we can obtain the corresponding diagram Fig|14b. Figll4b and Figll4ll from FigllSb. 
FigllSb and FigllSd. respectively. These diagrams are not drawn in Figll4[ They also contain collinear 
divergences. If the divergences survive in the end results, it implies that the factorization is broken. This 
needs to carefully be examined. 

We use k to denote the momentum carried by the gluon crossing the broken line. If the gluon is 
collinear, k has the patten: 



A;'^ ~ (1,A^A,A), A<1. 



(57) 



We use kg to denote the momentum carried by the gluon propagator with the short bar. The propagator 
has three terms in the light-cone gauge: 



n^kii n^'kf, 
.g^- + — ^ + — ^ 

77- * f\in Tl ' fXiri 



(58) 



In the above /i is the index contracted with that in the vertex left to the short bar, and v is contracted 
with that in the vertex right to the short bar. The first term will not give collinear divergence in Fig|14b. 
But, the second and third term will give collinear divergences with the collinear power-counting, because 
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the denominator of the terms is at order of A^, i.e., n ■ kq ^ X'^ derived from the on-shell condition 5{kg) 
with Eg. (1571). 

The propagator in Eq. ([58|) also appear in FigfTSb. The second term gives no contribution because of 
v{p)n ■ 7 = 0. The contributions from the first- and third term contain the colhnear divergences. It is 
easy to show that the divergence from the third term is canceled by that from the third term in Fig|14b. 
This also happens for other diagrams in FigJT3] in a similar way. Through explicit calculation we find 
that the divergence introduced by the second term in FigJ14k and FigJ14b are canceled by that in Fig|14b 
and Fig|14d. respectively. Therefore, only the collinear divergences in Fig|13l introduced by the first term 
in Eq. (l58p survive at the end, if we include all diagrams from FigJlSI and FigJ14l in the gauge n ■ G = 0. 
The diagrams in the light-cone gauge by changing the attachment of the collinear gluon in the right 
part of diagrams in Fig {13] do not contain collinear divergences. This has the implication for using the 
diagram expansion in the light-cone gauge, where one will have the uncanceled divergences from the cut 
gluon-propagator. With the method in Feynman gauge one will not have such divergences. 

Prom the above discussion the correct result is to obtain by taking only the first term in Eq. ()58p to 
calculate the diagrams in Fig|131 or by taking all in Eq. ()58p to calculate all diagrams inFiglTSland FigJ14[ 
We obtain: 
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With the results of relevant twist-3 matrix element in Eq.([5 
form: 



one can derive the following factorized 
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(60) 



For the (^f^-contributions there are also a SQP contribution, where one can obtain Tgib(x,0) from the 
^((/-contributions at one-loop. The SQP contribution in Wt is obtained by replacing hs with a gluon at 
one- loop. This contribution is in fact contained in the factorized from in Eq. (j56p . This is similar to the 
case in ^^-contributions for the SGP-contributions with Fig l8] discussed in Sect. 5.1. 

Combining all fiavors we obtain then the factorized SQP contributions as: 



Wt 



SQP vT^A'cgi Jx yi Jy y2 

5Q,+(6,6)E4/^?(y2m+(yi,o) + cSQ,_(6,6)Ee'/^;(y2)V(yi,o) 

[9] [q] 

+5q,,-+(6, 6) E elfq{y2)Tq+{-yi,0) + cSq,,-_(6, 6) E fq{y2)Tq-{-yi,0) 

[q] [q] 



(61) 



In comparison with the existing results in derived with the method of diagram expansion our results 
of SQP contributions are different. The difference is of an overall factor of —2. We note that the SQP 
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contribution is proportional to q^^ in the limit q'j_ — t- 0. Hence, it is not a leading contribution in the 
limit qI/Q"^ < 1. 

7. Summary 

We have studied the collinear factorization of SSA in Drell-Yan processes. To derive all perturbative 
coefficient functions at leading order of Og in the factorization, we have studied the scattering with multi- 
parton states, in which the helicity of the states are flipped. SSA in such a scattering is nonzero. This is 
in contrast to the scattering with a transversely polarized single quark. In this case SSA is always zero 
because of the helicity conservation of QCD for massless quarks. 

We have calculated SSA in the multi-parton scattering processes and the relevant twist-3 matrix 
elements of multi-parton states. By using the results from our calculation SSA has been factorized as 
convolutions of twist-3 matrix elements of the polarized hadron, parton distribution functions of the 
unpolarized hadron and perturbative coefficient functions. All perturbative coefficient functions of these 
contributions are derived here at the leading order of a^. In the factorization there are HP-, SGP- and 
SFP-contributions. From our results, we find that SSA at tree-level is factorized as the HP contributions. 
But the SGP- and SFP- contributions are from a class of one-loop contributions to SSA. These one-loop 
contributions contain collinear divergences and they can only be factorized with the soft-pole twist-3 
matrix elements in which one of the active patrons carries zero momentum. These soft-pole twist-3 
matrix elements are zero at tree-level but nonzero at one-loop. This results in that the perturbative 
coefficient functions of SGP- and SQP contributions are at the same order as those of HP contributions. 
Hence, in the collinear factorization there is a nontrivial order-mixing. Such an order-mixing does not 
happen in the factorization only involving twist-2 operators. 

It is interesting to note that at one-loop SSA contains divergences caused by exchanges of a Glauber 
gluon, as discussed in Sect. 5. The divergences are factorized with the soft-gluon-pole matrix elements. 
This is in contrast to the factorization of unpolarized cross-section only with twist-2 operators, where it 
is well known that the divergences from exchanges of Glauber gluons are canceled[38l[39ll30]. In the case 
of SSA studied here with twist-3 operators, such divergences are not canceled and need to be factorized. 
This will have some implications for the study of factorizations in the framework of soft collinear effective 
theories of QCDjU]. 

Our results for the collinear factorization of SSA in Drell-Yan processes agree with those derived 
with the method of diagram expansion, except the SQP contributions studied in Sect. 6. Comparing 
the method of the diagram expansion, we believe that it has advantages to use our method with multi- 
parton states for analyzing factorizations of SSA and for calculating higher order corrections, because the 
involved calculations are of standard scattering amplitudes. The approach we have taken here provides 
another way to derive the collinear factorization of SSA in various processes. It will be useful to solve 
the discrepancy between results for SSA in [42], where the momentum of a lepton in Drell-Yan processes 
is measured. It will also be useful for solving the discrepancy of evolutions of twist-3 matrix elements 
derived in [43 ^ 144 1 [45]. We leave these for future work. 

Note Added: During the preparation of the paper the results of the SGP-contributions with gluonic 
twist-3 matrix elements is reported in [46] . The results there agree with ours in Sect. 5.2.. 
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